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ABSTRACT 

We study the formation of an entangled state in a one-dimensional chain of three 
nuclear spins system which interact weakly through the Ising type of interaction and 
taking into account first and second neighbor interactions. We can get this entangled 
state using two pulses (it/ 2 and it pulses), and we study the efficiency of getting this 
entangled state as a function of the ratio of the second neighbor interaction coupling 
constant to the first neighbor interaction coupling constant (J'/ J). We found that 
for J' / J > 0.04, the entangled state is well defined. 
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1. Introduction 



The huge interest in quantum computation and quantum information was triggered 
by the polynomial time solution of the prime factorization problem (Shor's algo- 
rithm [1]), the fast data base searching (Graver's algorithm [2]), error correction 
codes [3], robust entanglement [4], and the teleportation phenomenum [5]. Almost 
any quantum system with at least two quantum levels may be used, in principle, 
for quantum computation. This one uses qubits (quantum bits) instead of bits to 
process information. A qubit is the superposition of any two levels of the system , 
called |0) and |1) states, * = C |0) + Ci\l) with |C | 2 + \Ci\ 2 = 1. The tensorial 
product of L-qubits makes up a register of lenght L, say \x) = ■ ■ ■ ,io), with 

ij = 0, 1, and a quantum computer with L-qubits works in a 2 L dimensional Hilbert 
space, where an element of this space is of the form ^ = J2C x \x), with ^ \ C X \ 2 = 1. 

Quantum computers of few qubits [6] have been in operations and have been used to 
explore quantum gates, entanglement, small number Shor's factorization, small data 
base Grover's searching, teleportation, error corrections, and cryptography How- 
ever, to make serious computer calculations one may requires a quantum computer 
with at least of 100-qubits registers, and we think that will be hopefully achived in 
a near future. One solid state quantum comuter model that has been explored for 
physical realization and which allows to make analytical studies is that one made 
up by one-dimensional chain of nuclear spins sytems [7] , where the Ising interac- 
tion among first neighbor spins allows to inplement ideally this type of computer 
with 1000-qubits or more [8]. One of the important phenomena we studied with 
this model was the entangled state formation [9]. In this paper, we consider second 
neighbor spin interaction in a chain of three nuclear spins system. We show that 
this allows us to implement an entangled state using two pulses (tt/2 and 7r), and 
we determine the threshold of the second neighbor interaction coupling parameter 
to get a well define entangled state. 



2. Equation of Motion 

Consider a one-dimensional chain of three equally spaced nuclear-spins system (spin 
one half) making an angle cos 9 = l/\/3 with respect the z-component of the mag- 
netic field (selected in this way to eliminate the dipole-dipole interaction between 
spins) and having an rf-magnetic field in the transversal plane. The magnetic field 
is given by 

B = (bcos(ujt + if), —bsm(ujt + <£>), B (z)) , (1) 

where b is the amplitude of the rf-field, B Q (z) is the amplitude of the z-component 
of the magnetic field, u and 92 are the angular frequency and phase of the rf-field. 
So, the Hamiltonian of the system is given by 

210 

h = - ]r /i k • B k - 2 jn j2 m+i - u'h j2 im + 2 , m 

k=0 k=0 k=0 

where represents the magnetic moment of the kth-nucleus which is given in terms 
of the nuclear spin as /U k = ^l{Ik^k->^k)-< being 7 the proton gyromagnetic ratio. 
B k represents the magnetic field at the location of the kth-spin (z = Zk). The 
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second term at the right side of (2) represents the first neighbor spin interaction, 
and the third term represents the second neighbor spin interaction. J and J' are 
the coupling constants for these interactions. This Hamiltonian can be written in 
the following way 

H = H + W , (3a) 



where Hq and W are given by 



U=o J 



and 



2 / .1 



(36) 



(3c) 



where u>k = lB (zk) is the Larmore frequency of the kth-spin, Q = 76 is the Rabi's 
frequency, and 1^ = 1% ± z/^ represents the ascend operator (+) or the descend 
operator (-). The Hamiltonian Hq is diagonal on the basis {|i2*i*'o)}> where ij = 0, 1 
(0 for the ground state and 1 for the exited state), 

#oK2*i*o) = E i2ilio |«2*l*'o) • ( 4a ) 
The eigenvalues Ei 2 i x i are given by 

E i2 ni = -^{(-ir 2 u; 2 +(-) l ^ 1 + (-iru;o+J[(-l)^ +ll + (-l) ll+t2 ] + (-l) lo+i2 J / } . 

(46) 

The term (3c) of the Hamiltonian (3a) allows to have a single spin transitions on the 
above eigenstates by choosing the proper resonant frequency, as shown in Figure 1. 
In this work, we are interested in the transitions |000) < — > |001) and 1 001) < — ► 1 101) 
which have the resonant frequencies 

to = loq + J + J' (5a) 

and 

uj = uj 2 + J - J' . (5b) 
To solve the Schrodinger equation 

iK-Qt = # * - (6) 
let us propose a solution of the form 

*(t) = £C fc (t)|fc> , (6) 

where we have used decimal notation for the eigenstates in (4a), Ho\k) = E/,\k). 
Substituting (6) in (5), multiplying for the bra (m\, and using the orthogonality 
relation (m\k) = 5 m k, we get the following equation for the coefficients 

7 

ihC m = E m C m + Y,Ck(m\W\k) m = 0,...,7. (7) 

fc=0 
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Now, using the following transformation 

C m (t) = D m {t)e- lE - t/h 



(8) 



the fast oscillation term E m C m of Eq. (7) is removed (this is equivalent to go- 
ing to the interaction representation), and the following equation is gotten for the 



coefficients D„ 



1 



iD m = ^ w mkD k e 



k=0 



where W m k denotes the matrix elements (m\W\k), and uj mk are defined as 



^mk 



E m — E k 

h 



(9a) 



(96) 



Eq. (9a) represents a set of sixteen real coupling ordinary differential equations 
which can be solved numerically, and where W mk are the elements of the matrix 
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(9c) 



where z is defined as z = e l M+v) ; an d z* is its complex conjugated. 



3. Numerical Simulations 

We start with the ground state, = 1 000), of the system and apply a 7r/2-pulse 
with tp = and with frequency to = ujq — J — J' to get the superposition state 

*i = ^(|000> + |001>) . (10) 

Then, we apply a 7r-pulse with 93 = and frequency uj = U2 + J — J 1 to get the 
entangled state 

* 2 = -L(|ooo)-|ioi>) . (11) 

The entangled state with plus sign can be gotten using a phase ip = ir. To solve 
numerically (9a), we select similar values for the parameters as reference 8 and 9. 
So, in units of 2tt x MHz, we set the following values 

w = 100 , wi = 200 , uj 2 = 400 , J = 5 , Q = 0.1 (12) 

The coupling constant J' is chosen with at least one order of magnitude less that 
J since in the chain of spins one expect that second neighbor contribution to be at 
least one order of magnitude weaker than first neighbor contribution, depending on 
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the interspace separation of the nuclei. In all our simulations the total probability, 
J2 \ Ck(t)\ 2 , is conserved equal to one within a precision of 1CT 6 . 

Figure 2 shows the behavior of Re Dq, Im Dq, Re -D5, and Im D§ during the two 
pulses (t = t = 7r/20 + tt/U) for the digital initial state and with J' = 0.2. We 
can see the formation of the superposition state after the first 7r/2-pulse and the 
formation of the entangled state (11) after the following 7r-pulse. Fig. 3 shows the 
behavior of the probabilities |Cfc| 2 , k = 0, . . . , 7 during the two pulses with the clear 
formation of the superposition state and the entangled state. Fig. 4 shows the 
behavior of the expected z-component of the spin, 

<^> = 5E(- 1 ) fc i c *(*)i 2 . ( 13fl ) 

fc=0 

(I!) = \{\Co\ 2 + |Ci| 2 - |C 2 | 2 - |C 3 | 2 + |C 4 | 2 + \C 5 \ 2 - |C 6 | 2 - |CV| 2 j , (136) 
and 

{I Z 2) = \jl\C k \ 2 -jz\Ck\ 2 , (13c) 

fc=0 fc=4 

during the two pulses. As one could expect = (if) = at the end of the two 
pulses due to the formation of the entangled state (11). The expected value of the 
spin is rotating in the plane (x, y) as is shown on Fig. 5. These transversal expected 
values are given by 



and 



(IS) = Re (C?C + C* 3 C 2 + C* 5 C 4 + C* 7 C 6 j , (I y ) = Im{. . .) , (14a) 
(If) = Re (c* 2 C + C 3 *d + C%C A + C 7 *C 5 ) , (if) = Im (...) , (146) 

(if) = Re (clCo + CtC x + CIC 2 + C 7 *C 3 ) , (i|) = Im (...) . (14c) 

To determine the range of values of J' (second neighbor coupling constant) for 
which the entangled state is well defined after the two pulses, that is, where the 
other resonances and non-resonant transition are canceled, we calclate the fidelity 
paramete [10] for this process, 

R — expected] ^numerical) j (1<^) 

where ^expected is our state (11), and ^ nume rical is the resulting wave function from 
our simulations. Fig. 6 shows the fidelity as a function of the ratio of the second 
neighbor interaction constant to first neighbor interaction constant, J'/ J. As one 
can see, for a value J'/ J > 0.04 one gets a very well defined entangled state. This 
means that the second neighbor interaction with even two orders of magnitud weaker 
than the first neighbor interaction, we can generate an entangled state in this system. 
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4. Conclusions and Comments 



We have studied the formation of an entangled state using two pulses in a chain 
of three nuclear spins system with first and second neighbor Ising spin interaction. 
The characteritics of the entangled state were determinated, and we found that the 
entanglement can be realized even for very weak second neighbor spin interaction 
(</' I J > 0.04). We consider that the coupling constant J' may play an important 
rolle in the so called 27r-method found in reference [8,9] to supress non-resonant 
transition in the chain of nuclear spins system because this parameter will enter in 
the detuning parameter ( A = (E p — E m )/h — uo). 
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Figure Captions 

Fig. 1 Energy levels and resonant frequencies of interest. 

Fig. 2 Entangled state formation, (1): Re D , (2): Ira D , (3): Re D 5 , (4): Ira D 5 
with J' = 0.2 

Fig. 3 Probabilities for J' = 0.2, (k): \C k {t)\ 2 for k = 0, . . . , 7. 
Fig. 4 Expected values (a): (b): {If), and (c): (If). 

Fig. 5 For J' = 0.2, expected values of the transversal components of the spin. 
Fig. 6 Real, Imaginary parts of the Fidelity, and its modulus. 



7 



References 

1. P.W. Shor, Proc. of the 35th Annual Symposium on the Foundation 

of the Computer Science, IEEE, Computer Society Press, N.Y. 1994, 124. 

2. L.K. Grover, Phys. Rev. Lett., 79 (1997) 627. 
L.K. Grover, Science, 280 (1998) 228. 

3. P.W. Shor, Phys. Rev. A 52 (1995) R2493. 

E.Knill, R. Laflamme, and W.H. Zurek, Science, 279 (1998) 342. 
A. Stean, Proc. R. Soc. London Se A, 452 (1996) 2551. 

4. F. Schmidt-Kaler, S. Guide, M. Riebe, T. Deuschle, A. Kreuter, G. Lancaster, 
C. Becher, J. Eschner, H. Haifner, and R. Blatt, 

J. Phys. B, 36 (2003) 623. 

C.F. Roos, G.P.T. Lancaster, M. Riebe, H. Haffner, W. Hansel, S. Guide, 
C. Becher, J. Eschner, F. Schmidt-Kaler, R. Blatt, 
Phys. Rev. Lett. 92, (2004) 220402. 

5. M. Riebe, H. HaiTner, C. F. Roos, W. Hansel, J. Benhelm, G. P. T. Lancaster, 
T.W. K6rber, C. Becher, F. Schmidt-Kaler, D. F. V. James, R. Blatt, 
Nature 429, (2004) 734. 

6. D. Boshi, S. Branca, F.D. Martini, L. Hardy, and S. popescu 
Phys. Rev. Lett., 80 (1998) 1121. 

C.H. Bennett and G. Brassard, Proc. IEEE international Conference on 

Computers, Systems, and Signal Processing, N.Y. (1984) 175. 

I.L. Chuang, N.Gershenfeld, M.G. Kubinec, and D.W. Lung 

Proc. R. Soc. London A, 454 (1998) 447. 

I.L. Chuang, N. Gershenfeld, and M.G. Kubinec 

Phys. Rev. Lett., 18 (1998) 3408. 

I.L. Chuang, L.M.K. Vandersypen, X.L. Zhou, D.W. Leung, and S. Lloyd, 
Nature, 393 (1998) 143. 

P.Domokos, J.M. Raimond, M. Brune, and S. Haroche, 
Phys. Rev. Lett., 52 (1995) 3554. 

J.Q. You, Y. Nakamura, F.Nori, Phys. Rev. Lett. ,91 (2002) 197902. 

7. Lloyd, Science, 261 (1993) 1569. 

G.P. Berman, G.D. Doolen, D.D. Holm, and V.I Tsifrinovich 
Phys. Lett. A, 1993 (1994) 444. 

8. G.P. Berman, G.D. Doolen, D.I. Kamenev, G.V. Lopez, and V.I. Tsifrinovich 
Phys. Rev. A, 6106 (2000) 2305. 

9. G.P. Berman, G.D. Doolen, G.V. Lopez, and V.I. Tsifrinovich 
quant-ph/9802015, quant-ph/9909032, Phys. Rev. A, 61 (2000) 062305. 

10 . A. Peres, Phys. Rev. A 30 (1984) 1610. 



8 



400 



TT- Hll> 



200 - 



N 

X 



(N 



- 



-200 



\t\f 



\f\f\f 



co 2 +J-J' 



I110> 



- I101> 



I100> 



I011> 



I010> 



I001> 



co +J+J' 



\f\t\f 



I000> 



-400 - 



Fig. 1 








10 



20 



30 



40 



50 



Fig. 2 




Fig. 4 



Ik MM 





